Bose Condensates in TOP Traps Exhibit Circulating Superfhiid Flows 
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For spin one atoms localized in a quadrapole magnetic field gradient, the atoms may be impeded from spin 
flipping their way out from the center of the trap by the application of a rotating uniform magnetic field. From 
a quantum mechanical viewpoint, such a trap for a Bose condensate is equivalent to having a superfluid in a 
rotating bucket. Vorticity is then expected to be induced in the condensate fluid flow without the application of 
any further external perturbations. 
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The two fluid modeljlj has been quite success- 
ful for describing liquid He 4 below the lambda 
temperature. The model asserts that the nor- 
mal fluid component can easily carry vorticity, 
but the superfluid component can undergo only 
potential flows with velocity v s = V<i>. However, 
it is presently understood that when superfluid 
if e 4 is located in a rotating container, vorticity 
can enter into the superfluid How in the form of 
"vortex lines" with a circulation k = § v s • dv 
quantized in units of kq = (2nh/M). 

In the light of recent progress in preparing 
Bose condensates in magnetic traps there 
has been considerable interest in whether or not 
vorticity can play an important role in the ex- 
perimental features of a mesoscopic superfluid. 
Our purpose is to point out that superfluid vortic- 
ity must already be present in those experiments 
which employ rapidly rotating magnetic fields, i.e. 
in the so-called "TOP trap" experiments. The de- 
tailed description of TOP traps will be reviewed 
in the work which follows. Here, we note that the 
implications of rotational vorticity for TOP trap 
experiments have not been previously explored. 

Consider (at first) the quantum dynamic be- 
havior of (say) superfluid He 4 in an arbitrarily 
shaped rotating bucket. Let H denote the Hamil- 
tonian of the superfluid in the bucket if the bucket 
were not rotating. The Hamiltonian H does not 
depend on time. However, since the bucket does 
rotate at angular velocity 

n = n(d9/dt), (I) 



where n is a unit vector along the axis of rotation, 
the fluid laboratory frame Hamiltonian actually 
develops a time dependence given by 

H{t) = S(t)HS*(t), (2) 

5(t)=eatp(-tn-Lfl(t)/n) > (3) 

where L is the total angular momentum of the 
fluid. All that is required to understand the above 
time dependent Hamiltonian H (t) for a fluid in a 
rotating bucket is that the angular momentum of 
the fluid is the generator of rotations]?]]. 

A crucial step is to eliminate the time depen- 
dence from the Hamiltonian H(t) by using the 
following quantum mechanical canonical transfor- 
mation (H(t) — ► Ti), 

H = S*(t)H(t)S(t) - ihS*(t)(8S(t)/dt), (4) 

to a frame rotating along with the bucket. Em- 
ploying Eqs.(l)-(4), one finds 

n = H-n-L. (5) 

Note that the effective Hamiltonian Ti in Eq.(5) 
does not depend on time if the angular velocity 
17 does not depend on time. We are thereby per- 
fectly justified^ in using Ti in Eq.(5) to describe 
a thermal equilibrium Bose fluid within a rotating 
bucket; The Gibbs canonical distribution is 

Pe9 =exp((F-H + n.-L)/k B l) (6) 

The notion of thermal equilibrium does not ap- 
pear so simple before the transformation from the 



laboratory frame (where H (t) depends on time) 
to the rotating frame (where TL does not depend 
on time). It is important to realize that Eq.(5) 
represents an entirely rigorous Hamiltonian which 
describes exactly why a rotating bucket at an- 
gular velocity l~i induces vorticity and net angu- 
lar momentum into the fluid contained within the 
bucket. 

The transition from a laboratory frame TOP 
trap time varying Hamiltonian H{t) to a rotating 
frame time independent Hamiltonian TL is a bit 
more subtle but will now be shown to be closely 
analogous. 

A TOP trap for triplet (spin one) atoms is 
simply a time dependent magnetic field B(r, t) 
constructed by superimposing a homogeneous 
rotating magnetic field B/j(i) together with a 
non-rotating inhomogeneous quadrapole mag- 
netic field Bg(r); In detail, 

B(r,t)=B h (t) + B Q (r), (7) 

where the quadrapole field is given in terms of the 
field gradient amplitude G, 

B Q (r)=G(r-3(n.r)n). (8) 

The homogeneous field B/j rotates about (and is 
normal to) the rotation axis unit vector n; In de- 
tail, n • Bo = and 

B h (t) = B cos(Qi) +n X B sin(fit). (9) 

A triplet state (spin one) atom interacts with the 
magnetic field B so that there exists three possi- 
ble (potential) energy levels, namely zero energy 
and ±7i7|B|, where 7 denotes the magnitude of 
the gyro-magnetic ratio. An atom in only one of 
these states can be trapped with the potential 

tf(r,t)=ft 7 |B(r,f)|, (10) 

The trap potential in cylindrical coordinates r = 
(p, ip, z) (where z — n • r) is determined by 

|B| = Gy/p 2 + b 2 + 2bpcos(<p - fit) + 4z 2 . (11) 

The magnetic length of the TOP trap is defined 
as b= \B \/G. 

For N atoms in a TOP trap, the time depen- 
dent Hamiltonian has the form 

H(t)= h ^)+ E v *- ( 12 ) 

l<j<Ar l<j<k<N 



The two body pair interaction on the right hand 
side of Eq.(12) is assumed to conserve total an- 
gular momentum, but the time dependent single 
particle Hamiltonian sum on the right hand side 
of Eq.(12) does not conserve total angular mo- 
mentum n • L; i.e. 

hi(t) = -(^)v 2 s +U(p j ,<p j -Slt,z j ) (13) 

has a potential which depends on the angle ipj, 
as well as pj and Zj. 

One may now solve the Schrodinger equation 
for the many body wave function x, 

in{ dX{ -2 , ''' ,t) )=m)x{..,r j ,...,tl (14) 
by looking for a solution of the form 
X(...,Pj,fj,Zj...,t) = *(..., p j ,ip j -£lt,z j ...,t){l$) 
yielding 

^•■■^■"■ f ) ) = ^(...,r jv ..,0. (16) 
The new Hamiltonian TL does not depend on time, 

H= E ~ h i+ E "i*. 

l<j<N l<j<k<N 

h j = -(^y j + uto) + i m( 1 L). (is) 

Eqs.(17) and (18) yield the Hamiltonian form 

H(TOP trap) =H-SI L = H- QL Z (19) 

which does not depend on time, does not conserve 
angular momentum but does represent the central 
result of this work. In the rotating frame of the 
TOP trap, the single atom potential does not de- 
pend on time 

U(r) = hjG^p 2 +b 2 + 2b p cos ip + 4z 2 , (20) 

and does not conserve the angular momentum 
component L z = n • L. 

In the previous theoretical literature concern- 
ing Bose condensates in TOP traps||(J, a time 
averaged Hamiltonian was employed 

( H )= E ( h i)+ E v ^ ( 21 ) 

l<j<N l<j<k<N 



where 
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tf + Ufazi) (22) 
with the time averaged potential 

U(p,z) = — U{p,ip-nt,z)dt. (23) 

2ti" Jo 

We note in passing that the time averaged TOP 
potential U(p, z) is not very well approximated 
by a simple anisotropic harmonic oscillator po- 
tential. More to our central point, there is no jus- 
tification for employing the time averaged Hamil- 
tonian for atoms in a TOP trap. Eqs.(17), 
(18) and (20) provide (in a mathematically rigor- 
ous fashion) the appropriate Hamiltonian Ti for a 
TOP trap, just as the general Eq.(5) has long pro- 
vided the appropriate Hamiltonian for superfiuid 
He 4 in a rotating bucket, u/ 
In the Gross-Piteavskiip 2 
gas model of a Bose condensate, the TOP trap 
induced order parameter of the condensate ought 
to obey the an equilibrium equation which follows 
from our above considerations; It is 
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i>(p, z ) = MV'(P) ¥>, z ), (24) 



where p is the chemical potential, and a is the 
two body scattering length. The ground state or- 
der parameter Eq.(24) is more than a little com- 
plicated. However, some of the implications of 
Eq.(24) may be stated with confidence. 

Only in the limit SI — > can the Bose conden- 
sate order parameter ip(r) be chosen to be real. 
For the TOP trap situation with fl ^ 0, the order 
parameter is complex 



ip(r) = V«s(r)exp (iAT$(r)/ft)) (25) 

yielding via Eqs.(24) and (25) a non- trivial super- 
fluid flow velocity v s = V<3> and a non-zero mean 
angular momentum. The velocity flow is deter- 
mined by a rotational quantum Bernoulli equa- 
tion 



p = Ml -|V$p - (fix r)-V$ 



(26) 
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If the TOP trap were not rotating, the small 
scattering length of two 87 Rb atoms (a ~ 50A), 
along with the small density of the condensate j| 
n s , would lead to a localization length within 
the trap potential p[ minimum determined by 
d ~ (h/M~/G) 1/3 ~ .15pm. This localization 
length d is much smaller than the magnetic length 
of b ~ 8 x 10 2 /iTO, i.e. d << b. If the conden- 
sate stayed within a distance d of the minimum of 
the potential U(r), then there would be frequent 
spin flips resulting in the loss of atoms from the 
trap. The rotating potential causes the conden- 
sate atoms climb a bit up the potential wall away 
from the center, acting in a manner closely analo- 
gous to particles in a centrifuge. Thus lifted from 
the potential minimum, the spin induced lifetimes 
of atoms in the trap are increased. 

To see how the mechanics of the rotational 
equilibrium for the condensate density n s may 
work, let us first make an analogy with a roulette 
wheel and a rolling metal ball found in many gam- 
bling casinos (not yet having completely entered 
the digital computer game simulation age). The 
mechanical roulette wheel turns and the ball is 
rolled into a grove of radius b. The ball rolls 
around this grove while the wheel spins. Along 
the grove, there are many minor potential min- 
ima, i.e. one minor minimum for every possible 
number that can be bet to win on the roulette 
wheel. The wheel continues to turn, and fi- 
nally the ball falls into one of the potential min- 
ima. The ball then rotates with the same angu- 
lar velocity as the roulette wheel achieving "rota- 
tional thermal equilibrium" at the winning num- 
ber. The point is that the rotational velocity of 
the ball relative to the rotational velocity of the 
roulette wheel cannot differ from zero for a very 
long time. Otherwise, the bets would never get 
settled. The ball does not in the final stages of the 
game feel the "Time average Of the Potential", 
i.e that time average in which the ball is contin- 
ually bouncing over the minor potential barriers 
during the initial stages of the game. The ball 
does finally settle down to a smooth rotational free 



energy minimum. 

In the TOP trap, the condensate is placed 
in a rotating potential grove of radius b. The 
Bose condensate can hardly stand still for very 
long while feeling only the "Time average Of the 
Potential" ; i.e. the condensate can only shake 
rapidly up and down the z-axis for a limited 
amount of time. Then the condensate then starts 
to flow smoothly around the grove. The critical 
angular velocity for forming one quantum of cir- 
culation, 

K Q = (2ir%/M) ~ 4.6 x 10~ 5 cm 2 /sec 

around to rotation axis cannot be very much 
larger than f2 c ~ (kq/o 2 ) which is (in order of 
magnitude) how slowly the minute hand of a 
clock rotates. This is a much lower rotational 
velocity then the TOP angular velocity[[l3) of 
(A/27r) ~ 7.5 x 10 3 Hz. One thus expects per- 
haps ~ 10 5 circulation quanta to flow around 
the axis of rotation. This circulation represents a 
very large number of quanta indeed. Unlike the 
roulette wheel there is but one potential minima 
per turn around the grove. The circulating Bose 
condensate, then flows up and down the poten- 
tial as it forms in its own little oval flow (as in a 
toroidal pipe) around this funnel potential. 

What happens to this flow of condensate 
around a toroidal pipe (so to speak) when the 
TOP trap is removed? The diameter of the pipe 
will increase and the tangential velocities to the 
flow around the ring will decrease (due to angular 
momentum conservation). The oval which is left, 
and displayed in an experimental picture could 
hardly be circular. It would be elliptic due to the 
displaced asymmetric axis of rotation. The above 
considerations are consistent with TOP trap ex- 
periments. 
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